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Abstract 

We study the conductance for the current through a finite interacting, pos- 
sibly disordered region attached to two noninteracting leads within the Kubo 
formalism applying the perturbation theory for the local Fermi liquid de- 
veloped by Yamada & Yosida and by Shiba to this finite region. Assuming 
the validity of the perturbation expansion in the Coulomb interaction and the 
time reversal symmetry for the normal scattering potential in the finite region, 
we find that the contributions of the vertex correction to the dc conductance 
disappear at T = if the currents are measured in the leads. Consequently, 
the conductance is expressed as g = (2e 2 /h) |t(0)| 2 . Here t(0) is the trans- 
mission coefficient for single-particle-like excitation at the Fermi energy and 
is introduced by using the Green's function. The results are generalized to a 
quasi-one-dimensional system with a number of scattering channels. 

PACS numbers: 72.10.-d, 72.10.Bg, 73.40.-c 
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I. INTRODUCTION 



The effect of the Coulomb interaction on the process of resonant tunneling through a 
small region has been a subject of considerable investigation in the past few years. For 
studying the transport through such small systems theoretically, it is sometimes necessary 
to formulate the expression of the conductance so that both the Coulomb interaction and 
quantum mechanical interference effects can be treated exactly, and several attempts to 
obtain a rather general formulation have been made so far. EH The purpose of this paper is 
to present one such attempt based on the Kubo formalism with the perturbation theory in 
the Coulomb interaction. 

We consider a system in which the Coulomb interaction and the scattering potential 
due to the disorder are switched on only for the electrons staying in a finite region of an 
infinite system, and apply the perturbation analysis used in the microscopic version of the 
Fermi liquid theory for the impurity Anderson modellHl to the present system. We find 
that the dc conductance for zero-temperature T = is written in a Landauer-type formBi 
if the assumptions shown in the following are satisfied. This result may be regarded as an 
extension of the relation shown in the noninteracting case by Economou & Soukoulis for one 
dimension and by Fisher & Lee for higher dimensions .1^0 The assumptions to be made are 

i) the Coulomb interaction and the normal spin-independent scattering potential are re- 
stricted for electrons in a finite spatial region, and two semi-infinite noninteracting 
leads are coupled to this region, 

it) perturbation expansion in the Coulomb interaction is valid, 

Hi) the single-particle-like excitation at the chemical potential fi does not decay, i.e., the 
imaginary part of the proper self-energy which is no longer diagonal in the wave-number 
indices is zero at T = 0. 

The assumption i) specifies the model, and means that the Hamiltonian does not contain 
a spin-flip-scattering due to the Kondo-type exchange interaction, the spin-orbit coupling, 



and so on. The noninteracting leads make the scattering problem of the present system 
well defined, and enable us to introduce the scattering coefficients making use of the asymp- 
totic form of the single-particle Green's function in the leads. Furthermore, the presence 
of the semi-infinite leads seems to make the Fermi-liquid description for low-energy states, 
the wave-length of which is much longer than the size of the central region L , probable. 
Under the assumption ii), we use the diagrammatic analysis for the perturbation in the 
Coulomb interaction considering the central region as an impurity with a finite volume,iHi 
and calculate the dc conductance by extracting the cu-linear imaginary part of the current- 
current correlation function following Shiba.i The assumption Hi) was originally introduced 
by Langer & Ambegaokar as a basic assumption in deriving the Friedel sum rule for interact- 
ing Fermi systems.B As Langer & Ambegaokar also mentioned, this assumption holds in the 
perturbation theory. It means that Hi) is not an independent assumption and deduced from 
the two assumptions i) and ii). We will show this explicitly by calculating the imaginary 
part of the self-energy following Yamada & Yosida.&S 

In this paper we consider mainly the system in one dimension. This is because the ex- 
tension to higher dimensions is possible following along almost the same lines. We show the 
outline of the extension in the last part of this paper. Owing to the presence of the non- 
interacting leads, the perturbation theory will probably make sense for low-energy states 
with long enough wave-length even in the one-dimensional case. The Tomonaga-Luttinger- 
liquid-like behavior seems to be valid for higher- energy region where the wave-length of a 
massless-boson-like excitation is shorter than the size of the central region Recently, 
the value of the renormalization factor K of the conductance in the Tomonaga-Luttinger 
liquid, g = Ke 2 /h (per spin), was discussed by several authorsEHl^l in relation to the ex- 
periment by Tarucha, Honda, and Saku.0 In the perturbation theory, the conductance can 
be expressed in a Landauer-type form, so that g becomes e 2 /h (per spin) when a resonant 
tunneling occurs, i.e., in the case the transmission coefficient for the single-particle-like ex- 
citation at the Fermi energy is unity. This result itself seems to be consistent with the 
result obtained by Shimizu with the phenomenological description of a Fermi liquid in one 



dimension. Ill 

We note that some of the ideas used in the present study were already applied to the 
current through a single Anderson- Wolff impurity00 in our previous work,0 where the 
expression of the conductance was studied for the use of numerical calculations with the 
quantum Monte Carlo method.0'@ The expression for the conductance through the single 
impurity has also been studied precisely with the Keldysh formalism by Hersfield, Davies 
and Willkins.il The conductance through a two-impurity Anderson- Wolff model has been 
formulated by Izumida, Sakai and Shimizu with a linear response theory for the use of the 
calculations with the numerical renormalization group method.@ 

In section ||, we introduce the transmission and reflection coefficients, and confirm that 
the coefficients for the single-particle-like excitation at the Fermi energy satisfy the unitarity 
condition. In section |HI| , we calculate the imaginary part of the proper self-energy with the 
perturbation theory in the Coulomb interaction, and show the assumption Hi) is deduced 
from i) and ii). In section [IV], we study the expression for the dc conductance g within the 
Kubo formalism, and show that the contributions of the vertex correction to g are zero at 
T = if the currents are measured in the noninteracting leads. Owing to this property, g 
can be expressed in terms of the transmission coefficient. We also discuss the role the vertex 
correction with relation to the back flow. In section |V|, we generalize the results obtained 
for one dimension in Sec. [I I 4IV| to the quasi-one-dimensional system which is finite along 
the transverse direction. 



II. SCATTERING COEFFICIENTS 

In this section, we introduce the transmission and reflection coefficients for a single- 
particle-like excitation, and show that the coefficients for the excitation at the Fermi energy 
satisfy the unitarity condition if the assumption Hi) holds. 

We consider an one- dimensional system which consists of three parts: a finite central 
region at < x < L, and two semi-infinite lead wires at — oo < x < and L < x < +oo . 
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The Hamiltonian is given by 
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where creates an electron with spin <r at position x. The spin- independent scattering 

potential V(x) and the Coulomb repulsion U(x, x') are restricted for the electrons staying in 
the central region; and U(x,x') = U(x',x). When the size of the central region L is of the 
order of the Fermi wave-length, the model is reduced to an Anderson-Wolff impuritylll'0 In 
this limit the perturbation theory is valid for all values of U. Our working hypothesis is that 
the ground state is changed continuously when the size L is increased from an atomic size 
to larger one owing to the presence of the noninteracting leads. Although it is not evident 
that the hypothesis is valid for whole region of the parameter space of the Hamiltonian, 
which may be valid at least for a finite region in the parameter space where U, V, and L are 
small enough. Unless, otherwise noted, we will be using units in which the Planck constant 
is h = 1. 

We now consider the transmission and reflection coefficients for the system described by 
the Hamiltonian Eq. ([J). For this purpose, we introduce the retarded Green's function 

G{x,x'-u + iQ + ) = -i J^dt ({^(x,t), ^o)})^ ^, (2) 

where 0(t) = e lHt Oe~ int , (■ ■ ■) denotes the thermal average Tr e~^ n ■ ■ ■ /Tie~ l3n with 
P being the inverse temperature 1/T, and the curly brackets denote the anticommutator. 
The spin index has been omitted from the left-hand side because the expectation value is 
assumed to be independent of whether spin is up or down. The Dyson equation for the 
Green's function is written, in the real space, as 

G(x, x'; z) = Gq(x, x'\ z) + / / dx\dx^ Gq(x, x\\ z) S(ari, £2; z) G(x2, x f ; z) , (3) 
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G (x,x';z)+ / / dx\dx 2 G (x, x\\ z) T(x\, x 2 ; z) G (x 2 , x'\ z) 
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(4) 



where Gq( 30 j OC - Z J IS the free Green's function corresponding to 7Yo, and the region of the 
integral is taken to be < xi, x 2 < L because the Coulomb interaction and the scattering 
potential are restricted for electrons in the central region. The proper self-energy S(xi, x 2 ; z) 
and the scattering matrix T(xi,x 2 ; z) contain both of the effects due to Tiu and TCy- Note 
that the functions introduced here are symmetric with respect to the permutations of x and 
x' because of the time reversal symmetry of 7i 

G(x,x';z) = G(x',x;z) 

T,(x,x';z) = Ti(x',x;z) ■ (5) 
T(x,x';z) = T(x',x;z) 

Therefore, the discontinuity of the Green's function at the real axis in the complex z-plane 
corresponds to the imaginary part 



G(x,x';e + iO + ) - G(x,x';e - i0 + )} /2i = Im G(x, x'\ e + i0 + ) 



(6) 



and the similar relation holds also for S(a?, x'\ z) and T(x, x';z). Specifically, the free retarded 
Green's function is obtained, for u > —p, as 



G (x,x'-u + iO + ) = -mp(u;)e iklx - x ' 1 , 



(?) 



where k = y2ra(u; + p), and p(oo) = m/irk is the density of states. 

When both x and x' belong to the leads, the asymptotic form of the full Green's function 
can be written explicitly, making use of Eqs. (f|) and (|7|): 



G(x,x';£ k + iO + ) = e lkx 1 - inpfa) T kk (£ k + i0+) {- mp ^ k )}e 



-i kx 



(8) 



for x' < and L < x; 



G(x, x'; 6 + *0+) = e ikx - inpfa) 7l fefc (£ fe + z0+) e^ kx { -mp{^ k ) } e 



-i kx' 



(9) 



for x' < x < . Here T k i k (z) is the Fourier transform of T(x, x'\ z) 



T k , k (u + iO + ) = [ L [ L dx 1 dx 2 e~ ik ' Xl T(xi, x 2 ; u + iO + ) e ikX2 , (10) 
Jo Jo 

and the frequency u is taken to be an on-shell value uj = ^ [= k 2 /2m — fi] in Eqs. (|8|) and 
(|). Physically, the retarded Green's function G(x, x'; + «0 + ) represents a propagation of 
a single-particle-like excitation; an additional electron for > 0, or a hole for — /j, < < 0. 
When £fc > 0, for instance, the right-hand side of Eq. (§) is relating to a probability amplitude 
for the process in which an electron created at x 1 < in the left lead is transmitted to the 
right lead with L < x. Similarly, Eq. (||) represents the process in which an electron created 
at x' < propagates as an incident wave e lkx or a reflected wave e~ lkx for x' < x < 0. 
Therefore it seems natural to define the transmission and reflection coefficients using Eqs. 
(1) and (|); 

t(&) EE 1 - Z7Tp(60 T fcfe (£ fc + *0+) 

r (£fc) = -inp(€k) T-kk(£k + i0 + ) 
With this definition, the unitarity condition |t(£fc)| 2 + l r (£fc)| 2 = 1 is expressed in the form 
of an optical theorem 

imr fcfc fe + zo + ) = [|r tt (6 b + io+)| 2 + |r_ w (& + io + )| 2 ] . (12) 

In the non-interacting system, the coefficients defined by Eq. flTT|) coincide with the scattering 
coefficients which are obtained from the asymptotic form of the wave function, and thus Eq. 
(fH2|) holds automatically. However, when the Coulomb interaction is switched on, it is not 
obvious whether Eq. (|T^) holds or not. We show in the following that the unitarity condition 
Eq. flT2"|) holds at T = for the single-particle-like excitation at Fermi energy = if the 
assumption Hi) mentioned in the Introduction holds for any positions x and x' in the central 
region as 

ImE(x,i'; i0 + ) = . (13) 

For simplicity, we temporarily write the integral equation Eq. (|3p in the operator form, 
G = Go + GoS G, and take the frequency for the operators to be z = i0 + . In this form the 
T-matrix is written as 



T = £ + £ G £ + £ G £ G £ + • • • . (14) 

Since Eq. QIBD means that £ is Hermitian at T = 0, i.e. £ = £*, the optical theorem can 
be derived following essentially the same way with that in the noninteracting case making 
use of an operator identity!! 

r - rt = (s + s g£ s + ■ • • ) (g - gJ) (£ + s g £ + ■ • • ) 

= Tt(G -Gj)T. (15) 
A diagonal element in the wave-number indices is written as 

-co 2^ Tfc '^ 0+) 27r ^ fe0 Tfc ' fe(Z0+) ' (16) 

where k is the Fermi wave number which satisfies = 0. Since the contribution of the 
integral in the right-hand side of Eq. ( |TED comes only from k! = ±k, the optical theorem 
Eq. (|T2p holds for the single-particle-like excitation at = 0. 

As already mentioned, the assumption Eq. (^) was introduced by Langer and Ambe- 
gaokar to derive the Friedel sum rule for interacting Fermi systems:@ it is expressed in the 
present case as 

I i-O c+oo 



An = / duo I dx Im 

7T J— oo J — oo 



G(x, X] uj + i0 + ) — Gq(x, x; oj + i0 + ) 



= -L Tr log [/+ (G -Gl)T] , (17) 

where the frequency for the operators in the second line is taken to be z = i0 + . We note 
that Eq. (|T^) can be derived with the perturbation theory in the Coulomb interaction and 
the calculation is performed in the next section. 

For the later convenience, we introduce here the Matsubara Green's function 

G(x,x'-ie n ) = - fdr (T t ^(x,t) ^ M (x',0)) e i£nT , (18) 

J 

where O m {t) = e TU Oe~ TH , and e n = (2n + l)7r//3. The Dyson equation for the Matsubara 
Green's function is obtained by setting z = ie n in Eqs. (|3|) and (|j). The free Matsubara 
Green's function is obtained as 
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Gq(x, x'\ is) = m — r — , (19) 

iK(ie) 



K,(ie) = v2m 



1/1 



(p(ie) = arg(/x + ie) . (21) 

where sgne is the sign function. The range of the value of phase ip(ie) is defined as — n < 
if < 7T. Thus in the limit £ — > ; ^(iO^) = for \i > 0, whereas y^O 1 * 1 ) = ±7r for \x < 0. 
Now we note some characteristic features in the limit e — > + . When /i > 0, 60(1, x'; i0 + ) 
is an oscillatory function of \x — a/|, and the singular contribution with the sign function in 
the right-hand side of Eq. ( |19D corresponds to the imaginary part of the retarded function 
for u = 0, i.e. lmG (x,x';iO + ), owing to the symmetry property Eqs. (|5|) and (§). On 
the other hand, when /1 < 0, Gq(x, x'; i0 + ) is a real and exponentially decreasing function 
of I a: — x'\, and the singular term disappears from the right-hand side of Eq. ( |T9| ) because 
exp[i(p(i0 ± )/2] sgnO 1 * 1 = i. 

III. IMAGINARY PART OF THE SELF-ENERGY 

In this section, we calculate the imaginary part of the proper self-energy up to the order e 2 
with the perturbation theory in the Coulomb interaction.!!!!! We show that the imaginary 
part is proportional to e 2 at T = 0, and thus Eq. (|I3D holds in the perturbation theory. 

As already mentioned, the self-energy satisfies the relation corresponding to Eq. (^), 

Im S(x, x'; e + i0 + ) = [ S(x, x'; e + i0 + ) - E(x, x'; e - i0 + ) ] /2i , (22) 

owing to the symmetry property Eq. (|5|). It means that the imaginary part of the retarded 
self-energy S(x, x'\ e + i0 + ) is related to the singular sgne n contributions in the Matsubara 
function S(x, x'\ ie n ) by the analytic continuation, and this correspondence holds for any x 
and x' . Therefore we can obtain the e 2 -term following the calculations in the case of the 
usual Fermi hquidlii in which the self-energy is a diagonal quantity. Yamada showed 
precisely that the singular contribution up to the order e 2 n comes only from diagrams which 



contain three intermediate Green's functions.Qc3 Fig. [1] shows the skeleton diagram for the 
simplest example, and the expression for the imaginary part is written after the analytic 
continuation as@ 

ImE (2) (i,i': e + iO + ) 



-°° de'de" 



, e'-e e' 

coth — — tanh — — - 

2T 2T 



, :" , e" + e-e' 

tanh — — — tanh 
2T 



oo (2tt) 

x /7^^' y') U{x, y) Im y; e" + z0+) 
xj 2 ImG(a;,x / ;e , + zO + )ImG(y,?/ / ;e + £"-£' + iO + ) 
- Im G(y, x'\ e' + i0 + ) Im G(x, y'\ e + e" - e' + i0 + ) j 



2T 



(23) 



Here G is the full Green's function which contains all the effects of the scattering potential 
and the Coulomb interaction. At T = 0, the hyperbolic functions in Eq. (^) are replaced by 
the sign functions as [sgn(e' — e) — sgn(e') ] [sgn(e") — sgn(£ // + e — e') ]. Thus the region 
of the integral for the frequencies are reduced to / de' I de", and the area of the integral 

JO Je'-e 

region becomes e 2 /2. Therefore the expression for small e is obtained setting the frequencies 
for the Green's functions in Eq. (|23|) to be zero; 

ImS (2) (i,x': e + i0 + ) 



-2 ,L r L 



2tt 2 Jo 



dydy U (x' , y') U (x, y) Im G(y', y; i0~ 



{2 lmG(x,x';iO + )ImG(y,y';iO + ) - lmG(y,x';iO + )lmG(x,y';tO + ) } 



(24) 



The imaginary part is proportional to e 2 for any x and x' because the e 2 -dependence is 
coming from the hyperbolic functions. It is convenient for proceeding next step to rewrite 
Eq. ( [ZIP by introducing the bare vertecies 7^ and 7^ [see Fig. |[ 



ImE (2) (i,i'; e + i0~ 



e 



1 ) 

2 ,L ,-L 



2tt 2 Jo Jo 



Jl dyi Im G(y 6 , y 5 ; i0 + ) Im G(y u y 2 \ i0 + ) Im G(y 3 , y 4 ] i0 + ) 



i=l 



1 



Tufa, Vs) V3, Vi) 7u(z/2, Va> a;') + - 7 TT (x, y 5 ; 2/3, yi) 7^(2/2, y±\ 2/6, a/) 



(25) 



where 
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7u(z4, x 3 ;x 2 ,x 1 ) = U(xi,x 2 )S(x 1 - x 4 ) S(x 2 - x 3 ) , (26) 
r r^(x^,x 3 ;x 2 ,x 1 ) = U(xi,x 2 )[5(x 1 - x 4 ) 5(x 2 - x 3 ) - 8(xi - x 3 ) 8(x 2 - x 4 ) ] . (27) 

We now include the higher order terms following Yamada,@'il which can be completed 
just by replacing the bare vertices 70-0-' with the full ones Y aa i [see Fig. |[ 

T aa , (x 4 , x 3 ; x 2 ,x 1 ) = r CTCT /. CT / CT (x 4 , x 3 ; x 2 , x 1 : 0, 0; 0, 0) . (28) 

This is because the singular e^sgne n contributions in the Matsubara function E(x, x';ie n ) 
can be obtained by summing up all the diagrams with the three intermediate Green's func- 
tions. Therefore, we can write down the expression which is exact up to £ 2 -term at T = 
by symmetrizing Eq. (^) with respect to the interchange of x and x' making use of Eq. 
(||), then adding the equivalent expression in which the set of variables for the integration 
{y5,y 3 ,Vi} and {y 6 , y 4 , y 2 } are interchanged, and finally replacing the bare vertices by the 
full ones; 



ImE(x,x'; e + i0 + ) 

6 

J[ dy { Im G(y 6 , y 5 ; i0 + ) Im G(y u y 2 ; i0 + ) Im G(y 3 , y 4 ; i0~ 



e 



2 r L pL 



lTT 2 JO JO i=1 

( x '> Ve'i 2/4, 2/2) T n {yi ,y 3 ;y 5 ,x) + T n (x, y 5 ; 2/3, 2/i) T n (y 2 , y 4 ; y 6 , x') 
+ r Ti( a:/ ' 2/5; 2/3, r n (y 2 , 2/4! ye,x) + T n (x,y 6 ;y 4 ,y 2 ) r n (y 1; y 3 ; y 5 , x') 
+ \ r n( x 'i 2/e! 2/4, 2/2) r n (yi, y 3 ; y 5 , x) + ^ T TT (x, y 5 ; y 3 , y x ) T n (y 2 , y 4 ; y 6 , x') 
+ 7} r TT ( x ' > V5', 2/3, yi) r n (2/2, 2/4; 2/6, + ^ r TT (x, y 6 ; y 4 , y 2 ) T v (y u y 3 ; y 5 , x') 



(29) 



Here we note that the full vertex function r CT4(J3;(T2(Tl (x 4 , x 3 ; x 2 , x\ : ie 4 ,ie 3 ;ie 2 ,iei) is anti- 
symmetric with respect to interchanges of the arguments (together with the spin suffixes) 
in the first or second pair: 1 and 2, or 3 and 4.il In the symmetrized expression, Eq. (p9"D, 
the eight terms in the bracket consist of four pairs each of which is connected by the sym- 
metry operation for the time reversal; (x 4 , x 3 ; x 2 , x\) ^ {xi ) x 2 ;x 3) x 4 ). Consequently, the 
imaginary part of the proper self-energy is proportional to e 2 and Eq. ( [13] ) holds at T = 
as long as the integral in Eq. (|29"D is finite, i.e., in the case the perturbation theory is valid. 
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IV. CONDUCTANCE 



In this section, we consider the expression of the dc conductance g applying the pertur- 
bation analysis! to the current-current correlation function. We show that the contributions 
of the vertex correction to g is zero at T = when the currents are measured in the nonin- 
teracting leads. The resultant expression, Eq. (|44l) , is written in terms of the transmission 
coefficient for the single-particle- like excitation at the Fermi energy, t(0), defined by Eq. 

When the electric field is applied to the central region, the dc conductance within the 
Kubo formalism is given by 

K(x,x';u; + iQ + )-K(x,x';iO + ) 
g = hm ; , 30 

where 

K(x,x';u + i0 + ) = i rdt([j(x,t),j(x',0)]) e^ +i0+ ^ , (31) 



fa) = Z 



t dip a (x) dipl(x) 

1>U*) — x ^- M*) 



(32) 



2mi 

As discussed by Fisher & Lee,0^ x and x' are arbitrary in Eq. ([$(]) and g is independent 
of the choice of the positions x and x' owing to the current conservation specific to one 
dimension. This feature can be seen explicitly in the Lehmann representation 

9 = Y~ E e~^(a\j{x)\a')(a'\j{x')\a)5{E a -E al ) , (33) 

act' 

where Z = Tre _/3W , |a) is the eigenstate of 7i, and E a is its eigenvalue. For obtain- 
ing this representation, we have used the relation caused by the time reversal symmetry: 
K(x,x';z) — K(x',x;z) and 

' K{x,x']u + i0 + ) - K{x,x';u-i0 + )] /2i = Im K(x, x'; u + i0 + ) . (34) 

In Eq. (pB]), the matrix element (a\j(x)\a') is independent of x because E a = E a >^ This 
can be confirmed by multiplying the equation of continuity dp(x,t)/dt +dj(x,t)/dx = by 
(ot| on the left and by \a') on the right, where p = ej^a^ti^a ■ 
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In the Kubo formalism, Eq. fl30|) , the dc conductance is obtained from the w-linear 
imaginary part of K(x,x';u) + i0 + ). In addition, Eq. ( |34] ) means that the imaginary part 
corresponds to the discontinuity of K(x,x'; z) at the real axis in the complex z-plane, i.e., 
the singular contribution of the form sgn(Imz). This type of singular contribution in the 
correlation function can be extracted using the perturbation analysis which was originally 
used by Shiba to derive the Korringa relation in Kondo alloys.il In the previous study, we 
applied this method to the conductance through a single Anderson- Wolff impurity.0 In order 
to carry out the analysis, we introduce the Matsubara function 

K(x,x';iis n ) = [ P dr{T T j M (x 7 r)j M (x',0)) e iv » T , (35) 
Jo 

where v n = 2im//3. The tu-linear imaginary part of K(x,x';uj + iO + ) is obtained from 
the | z/|-linear contribution of K(x,x';ii/) because \v\ = v sgnz/ is replaced by — iu by the 
analytic continuation iv —>■ u + iO + . The Matsubara function K(x, x'; iv) can be represented 
by a sum of the diagrams shown in Fig. |], and the corresponding expression for T = is 
written, replacing sums over discrete Matsubara frequencies by integrals, as 

K(x,x';iv) = K {a \x,x'-iv) + K {b) (x,x';iv) , (36) 

K^ a \x,x';iv) = - (— — ) {d/dx 4 - d/dxi}{d/dx 3 - d/dx 2 } 
\2miJ 

/+oo 
— G(x 3 ,x 1 ;ie)G(x A ,x 2 ;ie + iu) x , (37) 

K^ b \x,x']ip) = — (——] {d/dxi — d/dxA {d/dx 3 — d/dx 2 } 
\2miJ 

/+oo circle' r L 4 
77— /■•■/ Y[dyiG(y h xi;ie)G(x4,y4;ie + iv) 
rr.r -» (2vr) 2 Jo Jo f = \ 

xr fffrV v(2/4, 1/3; V2, Vi-i£ + iv, is'] it' + iv, it) 
xG(x3,2/3; ie') G(y 2 , x 2 ; ie' + iv) 



(38) 



In what follows we choose x to be in the right lead and x' to be in the left lead, and thus 
x'<0<yi<L<x with i — 1, . . . , 4. Then the asymptotic form of the Matsubara Green's 
functions in Eqs. ( |3"T| ) and fl3"8| ) can be written, making use of Eqs. (|j) and (|T^)-(PT|), as 

G(x,yf,ie) oc exp[z«(z£) (x — and G(x',yi]ie) oc exp[iK(ie) {yi — x')] . Therefore the 
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derivative with respect to x%, ... ,£4 can be done explicitly, and the expression valid for 
x' < and L < x is obtained as 

fT (a) (x,a;';iiz/) = - — ^/ — {k(ie + ii/) - k(«£:)} 2 G(x', x; ie) G(x, x'\ ie + iu) , (39) 
\2/mJ a J— 00 27r 

x + zv) — K,(ie) } G(y\, x\ ie) G(x, y^, ie + iv) 
xr^y^fe, Us, 2/2, yi : ie + ie'; ie' + iz/, ie) 

x{n(ie' + iu) — n(ie')} G(x' ,y 2 ;ie') G(y 3 ,x';ie' + iu) . (40) 



We can now extract the singular z/sgnz/ contribution from Eqs. (|39|) and (^0|). Shiba 
has shown precisely that this type of singular contribution comes only from the diagrams 
with two intermediate Green's functions which carry the same frequencies in the limit v — > 
0.1-0 The simplest example is the diagram Fig. ^ (a) in which the frequencies are assigned 
as G(ie)G(ie + iu). The z/sgnz/ contribution is originated from the sgne term in the 
Matsubara Green's function: for instance, see Eq. ( p0|) for the free Green's function. The 
full Green's function can be expressed, by separating the singular and regular parts, as 

G(x, x'\ ie) = G (x, x; ie) + %G (x, x'; ie) sgne . (41) 

Note that, when e — 0, G (x, x';0) and G '(x, x'; 0) correspond to the real and imaginary 
parts of the retarded Green's function G(x,x';i0 + ) owing to the symmetry property Eq. 
(§). The z/sgnz/ contribution of K(x,x';iu) due to the two intermediate Green's functions 
appears through the integration of the product sgn(e + u) sgn(e) over e ; 

/+00 r+00 
de A(ie) sgn(e + u) sgn(e) = / de A(ie) - 2 A(0) u sgn u . (42) 
-00 y — >o J— 00 

Here A(ie) is a smooth function which can be expressed in terms of G , G" , and so on, and 
the second term in the right-hand side appears as a result of the derivative d sgn(e + u)/du = 
2 5(e + u). Thus, the next thing we should do is to calculate the coefficient A(0) taking into 
account the vertex correction. In the perturbation theory, the diagrams which contain 
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the two intermediate Green's functions can be classified into four groups [see Fie. Ill In 
the figure, the two intermediate Green's functions which yield the singular contribution are 
marked with the cross, and the external frequency v is taken to be zero in the remaining part 
of each diagrams. One remarkable feature is that the contributions of the diagrams (bl)-(b3) 
become zero, i.e., the coefficient A(0) corresponding to these diagrams is zero. This is because 
the each of diagrams (bl)-(b3) contains at least one current vertex {n(ie f + iu) — K(ie')}, 
for which v can be simply taken to be zero. In contrast to (bl)— (b3), the contribution of the 
diagram (a) is finite because the two intermediate Green's functions and the current vertex 
carry the same the internal frequency e, and the product sgn(£ + v) sgnfe) appears also 
from the square of the current vertex [see Eq. (|39])]. Therefore, the finite contribution of the 
usgnu term comes only from the diagram (a), which can be calculated by substituting Eqs. 
(p0|), (pl|) and fl4ip into Eq. fl39|). Consequently, the dc conductance at T = is obtained, 
by reinserting fi, as 



e 2 1 



G{x,x';iO + ) , (43) 



nh [7rp(0) 

where x and x' is taken to be x' < and L < x as already noted. This expression can be 
rewritten in terms of the transmission coefficient t(0) defined by Eqs. (H) and (|TT|); 



9 = ^m I'- (44) 

These results Eqs. ([43]) and ([HP are generalized to a quasi-one-dimensional system with a 
number of channels in the next section. Furthermore, the present analysis can be applied to 
a tight-binding model and the results can also be generalized to the system on the latticed 
Before closing this section we consider briefly the reason why the contributions of the 
vertex correction disappear, and it seems to be understood in relation to the back flow in 
the leads. In the usual Fermi liquid theory the back flow can be represented in terms of the 
three point vertex function for the current ,0H which is one parts of Eq. ([38]) and can be 
written as [see Fig. H 

— / / dy 2 dy 3 T aa >. al(T (y 4 , y 3 ; y 2 , y x : is + iu, %e'\ is' + iu, is) 

-oo ZlT JO JO 

a 
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g 

x ^— {d/dx 3 - d/dx 2 } G(x 3 , y 3 ; is') G(y 2 , x 2 ; is' + iv) 



• (45) 

X2,Xs— >x' 



Specifically, the value A(x'; y±, y\, 0,0) for zero frequencies can be related to the renormal- 
ization of the current caused by the back flowBi Since we are considering the system 
without the translational invariance, a more precise definition of the back flow should be 
necessary for a general argument. Nevertheless, we can show the value A (a/; 2/4, y\\ 0, 0) to 
be zero when x' belongs to the noninteracting leads. It can be shown by performing the 
derivative with respect to it was done for obtaining Eq. (f40"D ; 

/+00 ds' fLrL 
— dy 2 dy 3 T a<jl . a , a {y^ y 3 ; y 2 , y 1 : is + iv, is'; is' + iv, is) 
-00 2,71 JO JO 

g 

x {niie + iv) — K{is')} G(y 2 , x ; is + iv) G(x', y 3 ; is') , (46) 

2m 

where ± is chosen to be + for x' < , and — for L < x'. Since the current vertex 
{n(ie' +iv) — n(is')} is zero for v = 0, the three point function becomes A(x'; y±, y±; 0, 0) = 0. 
In this sense, the current is not renormalized and the effect due to the back flow is absent in 
the noninteracting leads. This is caused by the fact that the asymptotic form of the Green's 
function in the leads is given by a simple superposition of the incident, transmitted, and 
reflected waves [see Eqs. (|j) and @]. 

V. EXTENSION TO HIGHER DIMENSIONS 

So far, we have considered the system in one dimension. In this section, we generalize 
the results obtained in the above to the quasi-one- dimensional system, which is finite in 
the direction perpendicular to the x-direction [the transverse direction p}. The extension 
can be done following along almost the same lines as those for the one-dimensional system. 
This is because the correspondence of the the imaginary part and the discontinuity of the 
Green's function holds also for the quasi-one-dimensional system, Eq. fl52|) , owing to the 
time reversal symmetry. So, the analysis of the Green's function used in Sec. 0-|IV| is also 
available here, and thus we will present only the outline of the extension. In Sec. |V A| , we 
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introduce the transmission and reflection coefficients for a number of scattering channels, 
Eq. fl55|), making use of the the single-particle Green's function, and examine the unitarity 
condition for the coefficients generalizing the discussion in Sec. [H|. In Sec. VB| , we consider 
the contributions of the vertex correction to the dc conductance in the quasi-one- dimensional 
system generalizing the calculations performed in Sec. [TV], and as a result the conductance 
is expressed in terms of the single-particle Green's function or the transmission coefficient; 
Eqs. © or 



A. Scattering Coefficients 

In this subsection, we generalize the definition of the transmission and reflection coeffi- 
cients discussed in Sec. || to the quasi-one- dimensional system. We now introduce a set of 
normal modes Xa{p) satisfying the Schrodinger equation for the transverse direction in the 
leads 

Xa(p) = £aXa(p) , (47) 

where a is the subband index, e a is the subband energy, and V c (p) is confinement potential 
for the transverse direction which is perpendicular to the x-direction. In what follows we 
choose the normal modes to be real, i.e., x a (.P) = Xl(p)- This is possible because of the 
time reversal symmetry of Eq. (f|7]), and is convenient for the discussion in the following. 
The free Green's function is diagonal in the subband indices 

G (r 1 r';z)=J2Xa(p)Xa(p')G^(x,x';z) , (48) 

a 

where r = (x, p), r' = (x', p'), and G^\x,x'; z) is the free Green's function for the 
one dimension, which is obtained by replacing the chemical potential fi by fi — e a in the 
right-hand side of Eqs. (0) or QTP|)-(PT[) . When the Coulomb interaction and the scattering 
potential are restricted for the electrons staying in the central region < x < L , the Dyson 
equation is written, by generalizing Eq. (||), as 



1 d 2 
-2^W + Vc(P) 
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G ba (x, x';z) = Gj, 0) (x, x'\ z) 5 ba 



L r L 



+ 11 dx 1 dx 2 G { b '(x,x 1 ;z)T ba (x u x 2 ;z)G^\x2,x';z). (49) 
Jo Jo 

Here the subscript ba denotes the subband indices defined by 

G ba (x, x'; z) = jjdpdp' Xb(p) G(r, r'; z) Xa(p') ■ (50) 

Because of the time-reversal symmetry, the relation corresponding to Eq. (|5|) holds also for 
the quasi-one-dimensional system as G(r', r; z) = G(r, r'\ z). Thus using Eq. (|50~D , we obtain 



G ab (x,x';z) = G ba (x',x;z) 



(51) 



Therefore, as in the one- dimensional case, the discontinuity of the Green's function at the 
real axis in the complex z-plane corresponds to the imaginary part; 



G ab (x, x'\ e + z0 + ) — G ab (x, x'\ e — i0 + ) /2i = lmG ab (x, x'\ e + iO 



(52) 



Owing to this property, we can obtain the imaginary part using the perturbation analysis 
for extracting the singular contributions. 

When both x and x' belong the leads, the asymptotic form of the full Green's function 
can be written explicitly making use of Eq. ([49]). So, the equation corresponding to Eq. @ 
is written by taking the values of x' and x to be x' < and L < x ; 



G ba (x,x';i0 + ) 



D i ktx 



5ba %a:k b k a {i0 

V b 



v a 



(53) 



where the frequency is chosen to be z = i0 + , k a = y2m( / u — e a ) is the wave number of 
single-particle-like excitation at Fermi energy, and v a = k a /m is the velocity. In Eq. (|53|) 
the initial and final states, a and b, are assumed to be propagating modes with real values 
of the wave number. T ba ± b k a (z) is the Fourier transform of % a {%i,%2'-i z) with respect to x\ 
and X2 [see Eq. fllTj)]. Similarly, the equation corresponding to Eq. (||) is written by taking 
the positions to be x' < x < ; 



G ba (x,x';i0~ 



5 ba e ikaX 



— %a-.-k b k a {i0 + ) e 

Vb 



-i ki,x 



I 

V a 



(54) 



Since the retarded Green's function G& a (x, x'\ i0 + ) represents a propagation of a single- 
particle-like excitation at the Fermi energy as discussed in Sec. [H], it seems natural to define 
the scattering coefficients making use of Eqs. ([53]) and fl5~4|) as 



% 

Vb 



$ba — — %a:k b k a {iO + ) 



n a = -— T ba: _ fei)fca (20 + ) J — , (56) 
v b V U " 

where a factor \Jvb/v a is introduced for the normalization. Then the unitarity condition, 
J2b [ \ ha\ 2 + \ r ba\ 2 ] = 1, is expressed in the form 

lmT aa:kaka ^0 + ) = -E'^[lW a (^ + )r + |T fea: _ fc6fea (^0 + )| 2 ] , (57) 

where the primed sum runs over the propagating modes. At T = 0, the unitarity condition 
holds for the excitation at the Fermi energy if the imaginary part of the proper self-energy 



lmE ab (x,x'; i0 + ) = . (58) 

As in the one-dimensional case, Eq. ( p7|) can be derived by using the operator identity Eq. 
(|15|). In the present case the matrix element of the operator T is defined by (b, kb\T\a, k a ) = 
%a-.k b k a (i0 + ) , and the product of the operators is understood to be an integral for the motion 
in the x direction and a sum over the subband index. The operator for the proper self-energy 
S is Hermitian owing to Eq. (|58|). Thus Eq. (|57| ) can be obtained from the diagonal element 

of the identity in the wave-number indices, which is written by generalizing Eq. fll~6|) as 

/+oo Ah' 

(59) 

Because of the energy conservation due to the delta function, the sum over b is restricted 
to the propagating modes and thus we obtain Eq. (|57|). We note that Eq. ( |58"D can also 
be derived in the perturbation theory by extending the calculation of the imaginary part 
performed in the Sec. |TJ; by including sums over subband indices for intermediate states 
and using Eq. (|52|) . 
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B. Conductance 



Next we consider the expression of the conductance by extending the calculation per- 
formed in Sec. [TV[ In the present case, the operator for the current in the x direction is 
given by 



J(x) = £ Jdp 



2mi 



(60) 



dx dx 

When the electric field is applied to the central region along the x direction, the dc conduc- 
tance within the Kubo formalism is given by 



K(x, x': uo + i0+) - K(x, x': i0 + ) 
g = hm : , 

K(x, x'; u + i0 + ) = i dt([ J(x, t) , J(x', 0) ] ) e l ^ +i0 + > * . 
Jo 



(61) 
(62) 



where x and x' are arbitrary owing to the current conservation. In the Kubo formalism 
the dc conductance is obtained from the cu-linear imaginary part of K(x,x';uj + i0 + ). In 
addition, the cu-linear imaginary part corresponds to the v sgn v term in the Matsubara 
function K(x, x'; iv) because the relation Eq. (|3~4] ) also holds for K owing to the symmetry 
property K(x, x'; z) = K(x', x; z). Therefore, the conductance is obtained by extracting the 
vsgnv term using the diagrammatic analysis. At T = 0, the equations corresponding to 
Eqs. (|36|)-(|38|) are written as 

(63) 



K(x,x';iv) = K {a) (x,x';iv) + K (b \x,x';iv) 



K {a) { 



( e 



\2mi 



{<9/<9x4 — d/dxi} {d/dxs — 8/8x2} 

/+00 
— G ab {x 3 , xi, is) G ba {x 4 , x 2 ; is + iv) 
-00 Ztt 



o ab 



(64) 



K^ b \x,x';iv) = — ( — ^ J {8/8x A — 8/dxi}{8/dx 3 — d/dx 2 } 

\ / I Lb / 



/+00 dsdE* rL rL ^ 
77TT2 /'"/ H d Vi G aib(yi,xi;is) G hai (x A ,y A ;ie + 

mr ao ^ ~°° K Zlr ) J ° J ° i=l 

0304 

xTatvT 1 (2/4, V3] V2, Vi-i£ + iv, is'; is' + iv, is) 
xG aa , i (x 3 ,y 3 ;is')G a2a (y2,x 2 ;i6' + iv) xi>x ^ x ■ 



IV) 



(65) 
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Here superscripts of the vertex function denote the subband indices defined by 

r 4 

(66) 

In what follows we choose x to be in the right lead and x' to be in the left lead, and thus 
x f <0<yi<L<x with i = 1, . . . , 4. Then, the derivative with respect to X\ . . . , x 4 can 
be performed explicitly making use of the asymptotic form of the Green's function, and the 
equations corresponding to Eqs. (|39| ) and ( f4"DD can be obtained as 

/ g \ 2 r+oo ds 

K (a >(x,x';iv) = - [ — ) / — {Kb(ie + iu) - K b (is)} {K a (ie + iu) - K a {ie)} 

\2mJ a ab J -co 27T 

xG a b(x' , x;is) Gba{x,x';ie + iv) , (67) 
2 r+™dede' r L r L A 



arr' „h ^ - X (2n) 2 Jo Jo , , 

a 3 a 4 

x{n b (is + %v) - K b (ie)}G aib (yi, x; ie) G ba4 {x, y 4 ; is + iv) 
xT a a 4 a ^ ai (y 4 , y 3 ; y 2 , y x : ie + iv, ie'; ie' + iv, ie) 

x{n a (ie' + iv) - K a (ie')} G aa3 (x' ,y 3 ;ie') G a2 a(y 2 ,x']ie' + iv) . (68) 
Here K a (ie) denotes the quantity which is given by replacing /i by /x — e a in Eqs. (p0|) and 



At T = 0, the contributions of the vertex correction, Eq. (|68|), to the dc conductance 
becomes zero as that in the one-dimensional case. This is because the current vertex {K a (ie'+ 
iv) — K, a (ie')} is zero for v = 0. This feature can also be understood in terms of the three 
point vertex function for J(x'), which represents the renormalization of the current or the 
effects of the back flow as discussed in the last part of Sec. [TV] [see Fig. ^jfl. Consequently, 
the dc conductance is determined by the v sgn v term in the right-hand side of Eq. (|67f) , and 
the expression of the dc conductance for T = is obtained by taking the asymptotic limit 
x' — > — oo and x — > oo in order to suppress the contribution of the evanescent modes with 
imaginary values of the wave number,0 and reinserting K\ 

g = —r J2v a v b G ba (x,x';i0 + ) , (69) 

nfl ba 
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where the primed sum runs over the propagating modes. This expression can be written in 
terms of the transmission coefficient t ab defined by Eqs. ( |53| ) and (^5j) as 




ba 



These are the results obtained with the perturbation theory. 

VI. SUMMARY 

We have studied the conductance for the current through a finite central region which 
is connected to semi-infinite noninteracting leads on the left and the right. Assuming the 
Coulomb interaction U and the normal scattering potential with time reversal symmetry V 
can be switched on only for the electrons in the central region, we have stressed the study of 
the relation between the conductance in the Kubo formalism and the transmission coefficient 
introduced by making use of the single-particle Green's function. In order to calculate the 
contributions of the vertex correction for the dc conductance, we have used the diagrammatic 
analysis for the perturbation expansion with respect to the Coulomb interaction, which was 
originally applied to the impurity Anderson model for the Kondo alloy. At T = 0, the 
contributions of the vertex correction to the dc conductance become zero if the currents are 
measured in the leads. This is caused by the fact that the asymptotic form of the Green's 
function in the leads is given by a simple superposition of the incident, transmitted, and 
reflected waves. Consequently, the conductance is expressed in a Landauer-type form using 
the transmission coefficient for single-particle-like excitation at the Fermi energy [see Eq. 
(fPl)l. At T = 0, the scattering coefficients are confirmed to satisfy the unitarity condition 
owing to a Fermi-liquid property Im £(x, x'; i0 + ) = 0, and this property can be derived with 
the perturbation theory. The results are generalized to a quasi-one-dimensional system with 
a number of scattering channels [see Eq. (f70|)1. The analysis used in the present study can 
also be applied to a tight-binding model, and the similar results can be derived for the model 
on a lattice S 
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Throughout the present study, we have assumed that the perturbation expansion in the 
Coulomb interaction is valid. For instance, in the limit the size of the central region L — > 0, 
the model is reduced to an Anderson- Wolff impurity and the perturbation theory is valid 
for all values of U. Although it is not evident that the perturbation theory is valid for all 
values of U, V, and L in the general case, the hypothesis may be valid at least for a finite 
region of the parameter space of the Hamiltonian where U, V, and L are small enough. 
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FIGURES 



FIG. 1. Second order skeleton diagrams for the proper self-energy x'; ie n ). 
FIG. 2. Second order skeleton diagrams in terms of the bare vertex function 7<j CT '(£4, X3; X2, xi). 
FIG. 3. Full vertex function T a4aa - CT2CT1 (x4, X3; X2, x\ : 1E4, ie^; 1E2, i£i) ■ 

FIG. 4. Current-current correlation function K(x, x'; iu). 

FIG. 5. Diagrams for the z/sgn v contribution in K \x , x' '; iu) . The two intermediate Green's 
functions which yield the singular contribution are marked with the cross. 

FIG. 6. Three point vertex function for the current A(x'; 2/4, y\\ is + iu, is). 



26 




Fig. 1, A. Oguri 




Fig. 2, A. Oguri 




Fig. 3, A. Oguri 
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Fig. 4, A. Oguri 




Fig. 5, A. Oguri 




Fig. 6, A. Oguri 



